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ABSTRACT 

A new control design methodology is introduced for multi-input/multi-output systems with 
unstable open loop plants and saturating actuators. A control system is designed using well known 
linear control theory techniques and then a reference prefilter is introduced so that when the 
references are sufficiently small, the control system operates linearly as designed. For signals large 
enough to cause saturations, the control law is modified in such a way to ensure stability and to 
preserve, to the extent possible, the behavior of the linear control design. 

Key benefits of this methodology are: the modified feedback system never produces 
saturatin • control signals, integrators and/or slow dynamics in the compensator never windup, the 
directional properties of the controls are maintained, and the closed loop system has certain 
guaranteed stability properties. 

The advantages of the new design methodology are illustrated in the simulation of an 
approximation of the AFTI-16 (Advanced Fighter Technology Integration) aircraft multivariable 
longitudinal dynamics. 
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1, Introduction 

In general the presence of saturations introduce stability and performance problem to closed 
loop multivarible control systems. One of those problems is the fact that control saturations alter 
the direction of the control vector. Each saturation element operates on its input signal 
independently of the other saturation elements; as it has been shown in the performance analysis in 
[10] and [11] this can disturb the direction of the applied control vector. Consequently, erroneous 
controls can occur causing degradation of the performance of the closed loop system over and 
above the expected fact that output transients will be "slower". 

Another performance degradation occurs when a linear compensator with integrators is used 
in a closed loop system and the phenomenon of reset-windup appears. During the time of 
saturation of the actuators, the error is continuously integrated even though the controls are not 
what they should be. The integrator, and other slow compensator states, attain values that lead to 
larger controls than the saturation limits. This leads to the phenomenon known as reset-windup, 
resulting in serious deterioration of the performance (large overshoots and large settling times.). 

The effects of the saturations on the closed loop stability of the control system are well 
known . When the open loop plant is unstable one can only guarantee local closed loop 'stability 
and there are references and disturbances such that when applied the closed loop system becomes 
unstable. 

Many attempts have been made to address these problems for SISO systems, but a general 
design process has not been formalized. One way to design controllers for systems with bounded 
controls, would be to solve an optimal control problem; for example, the time optimal control 
problem or the minimum energy problem etc. The solution to such problems usually leads to a 
bang-bang feedback controller [1]. Even though the problem has been solved completely in 
principle, the solution to even the simplest systems requires good modelling, is difficult to calculate 
open loop solutions, or the resulting switching surfaces are complicated to work with. For these 
reasons, in most applications the optimal control solution is not used. 

Because of the problems with optimal control results, other design techniques have been 
attempted. Most of them are based on solving the Lyapunov equation and getting a feedback which 
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will guarantee global stability when possible or local stability otherwise [2]-[3]. The problem with 
these techniques is that the solutions tend to be unnecessarily conservative and consequently the 
performance of the closed loop system may suffer. 

Attempts to solve the reset windup problems when integrators are present in the forward 
loop, have been made for SISO systems [4]-[6]. Most of these attempts lead to controllers with 
substantially improved performance but not well understood stability properties. As part of this 
research, an initial investigation was made on the effects on performance of the reset windups for 
MIMO systems [8] showing potential for improving the performance of the system. A simple case 
study was also recently conducted on the effects of saturations to MIMO systems where potential 
for improvement in the performance was demonstrated [9]. 

A new design methodology has been introduced in [11] for designing MIMO control systems 
for stable open loop plants with multiple saturations. Here a systematic methodology is 
introduced to design control systems with multiple saturations for unstable open loop plants. The 
idea is to design a linear control system ignoring the saturations and when necessary to modify that 
linear control law. When the exogenous signals are small, and they do not cause saturations, the 
system operates linearly as designed. When the signals are large enough to cause saturations, the 
references are modified in such a way to preserve ("mimic") to the extent possible the responses of 
the linear design. Our modification to the linear compensator is introduced at the error via a 
Reference Governor (RG). 

2. Performance Analysis 

Without loss of generality one can assume that each element u^t) of the control vector u(t) = [ 
T 

Uj(t) . . . tip(t)] has saturation limits ±1 and the saturation operator can be defined as follows: 


sat(u.(t)) 


1 

U(t) -1 <Uj(t) < 1 
-1 u.(t)<-l 


( 2 . 1 ) 
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Figure 2.1 shows the closed loop system with the saturation element f(u(t)) (f(u(t)) = 
[sat(uj(t) . . ,sat(u m (t)] T ) at the controls. The compensator K(s) is designed using linear control 
system techniques and it is assumed that the closed loop system without the saturations (the linear 
system) is stable with "good" properties. 


dj(t) d Q (t) 



Figure 2. 1 : The closed loop system 


There are well developed methods for defining performance criteria and for designing linear 
closed loop systems which meet the performance requirements. It would then be desirable, 
whenever the closed loop system operates in the linear region, to meet the a priori performance 
constraints (because it easy to define them and easy to design control systems satisfying these 
constraints). When the system operates in the nonlinear region new performance criteria have to be 
defined and new ways of achieving the desired performance must be developed. 

There are two major problems that multiple saturations can introduce to the performance of 
the system: (a) the reset windup problem, and (b) the fact that multiple saturations change the 
direction of the controls. 

When the linear compensator contains integrators and/or slow dynamics reset windups can 
occur. It is obvious that if the states of the compensator were such that the controls would never 
saturate, then reset windups would never appear. More information about the performance 
degradation of the feedback system caused by the saturations is given in [12]. 

To solve the performance problem a nonzero operator O 2 is applied in the reference 


signals and it will be called Reference Governor (RG). 
e = 0 2 r - y 


(2.2) 
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The nonzero operator is chosen, if possible , so that the control u(t) never saturates, i.e. 
llu^Hoo < 1, for any reference and/or disturbances. It is also desired that the 0 2 r signals are 
"close" to the r(t) signals. In section 4 the RG operator will be defined in detail. 


3. Mathematical preliminaries 

This section is an introduction to the new design methodology. Some necessary mathematical 
preliminaries will be given and a basic problem will be introduced. The basic problem will be 


solved and it's solution will lead to the new control design methodology. 

Consider the following linear time invariant system 

x(t) = Ax(t), x(0) = x 0 A e R nxn , x(t) e R n (3.1) 

y(t) = Cx(t) C e R mXn , y(t) e R m (3.2) 

y(x 0 ,t) = Ce At x 0 (3-3) 

where e At is the state transition matrix (matrix exponential) for A 

Definition 3.1 : The scalar-valued function g(x) is defined as follows: 

g(x 0 ): R n -* R, g(x 0 ) = ly(x 0 ,t)ll oo (3-4) 

Definition 3.2 : The set P g is defined as: 

P g = { [x,v] : xe R n , ve R, v > g(x) } (3.5) 

Definition 3.3 : Ba,c is the set of all xe R n with 0 < g(x) < 1, i.e. 

B A ,c = (x : 0 < g(x) <1} (3.6) 


Suppose that the system (3.1)-(3.3) has an initial condition x 0 e B AjC . From this definition 
we see that for such an initial condition the output of the system, y(t), will satisfy lly(t)lloo ^ 1. For 
neutrally stable systems the function g(x), the set P g and the set B Ai c have the following 
properties. 

(a) The set P g is a convex cone. 

(b) The B \,c set is symmetric with respect to the origin and convex. 

(c) .The function g(x) is continuous and even. 

(d) The function g(x) is not necessarily differentiable at all points in R n 
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Definition 3.4 [7] : The upper right Dini derivative is defined as 

+ f(t)-f(t 0 ) 

D f(t Q ) = lim sup — — 


The D + f(to) is finite at tg if the function f satisfies the Lipschitz condition locally around to [7], 
Note that the function g(x) defined by Definition 3.1 satisfies the Lipschitz condition locally if and 
only if the system (3.1)-(3.3) is neutrally stable. 

Theorem 3.1 171 : Suppose that f(t) is continuous on (a,b), then f(t) is nonincreasing on (a,b) 
iff D + f(t)<0 for every te (a,b). 



Figure 3.1: Visualization of the function g(x) and the sets P g and Ba,c- 



Assume that a stable linear system is defined by the following equations 

x(t) = Ax(t) + Br(t) A e R nxn , B e R nxm (3.8) 

y(t) = Cx(t) Ce R mxn (3.9) 

The goal here, is to keep the outputs of the linear system bounded (i.e. Iy ; (t)l < 1 V t, i) for 
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any r(t). In reference [1 1] a time-varying gain was introduced and this problem was solved 
completely. Here a time-varying-rate will be introduced and a different solution will be obtained. 
One can modify the inputs to the system r(t) to r^(t) with a time-varying rate operator, such that for 
any input r(t) the system output y(t) remains bounded. The new system can be defined as follows 
(also shown in figure 3.2). 


x(t) = Ax(t)+Br^(t) 

(3.10) 

^(0 = M-(t) ( r(t) - r^(t) ) 

(3.11) 

y(t) = Cx(t) 

(3.12) 


r R (t) 

I 


..A... 



The Basic Problem: 

At time t 0 find, if possible, the maximum time-varying rate p.(t 0 ), 0 < (i(t 0 ) 
<<», such that Vr(t), t > t 0 3 |i(t), t > t 0 such that the output will satisfy |y s (t) j < 1 


V i, t > t n 


Define the following auxiliary system 
x(t) = Ax(t)+Br^(t) 
r^(t)=|i(t)e r (t) 

i y(t) = Cx(t) 

and with x a (t) = [ x(t) r^(t) ] T one can obtain the augmented system 
x a (t) = A a x a (t)+B a p(t)e r (t) A a e R 


B a e R 


(3.13) 

(3.14) 

(3.15) 


n+m x n+tn 


n + m x m 


(3.16) 
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y(t) = C a x a (t) C a eR rnxn+m (3.17) 

To obtain the solution to the basic problem we define a function g(x) for the system (3.16)- 
(3.17) for B a = 0 and a set Ba,c as described previously. 

g(x a0 ): R n+m ^R, g(x 0 ) = Iy(x a0 ,t)H oo (3.18) 

x a (0) = x a0 , (3.19) 

B A ,c = {x: g(x) < 1 } (3.20) 

The function g(x) is finite for all xe R n since it is assumed that the system (3.8)-(3.9) is 
stable and consequently the system (3.16)-(3.17) is neutrally stable. 

Construction of fit): 

For every time t choose |i(t) as follows 

a) if x a (t)e IntB Ai c then |l(t) = <*> which implies that r(t) = r^(t) (3.21) 

b) if x a (t)e BdB Ai c then choose the largest |i(t) such that (3.22) 

0 < pt(t) < oo 


g(x a (t)+e[A a x (t)+B a fi(t)e (t)]) - g(x (t)) 

lim sup < 0 

e->o e 


(3.23) 


or for the points where g(x) is differentiable choose the largest |i(t) such that 
0<^t(t)<oo (3.24) 

Dg(x a (t)) [ A a x a (t)+B a |i(t)e r (t) ] <0 V t > 0 (3.25) 

where Dg(x a (t)) is the Jacobian matrix of g(x a (t)) as in definition 3.2. 

c) if x a (t)g B Ai c then choose p(t), 0 < )i(t) < «> such that the expression (3.23) is 
minimum. 

Theorem 3.2 f 1 21 : For the system given in (3. 16)-(3. 17) and the given construction of p(t) 
the following is always true Vx a e R n+m . 


g(x (t)+e[A x (t)]) - g(x (t)) 

lim sup ^ — <0 

£— e 


(3.26) 


and at the points where g(x) is differentiable 

Dg(x a ) A a x a < 0 Vx a e R n+m (3.27) 

where Dg(x a (t)) is the Jacobian matrix of g(x a (t)). 

The construction of |i(t) is always possible because of theorem 3.2 , namely one can choose 
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p(t) = 0 Vt and the inequality (3.23) is always true. 

Lemma 3.1 f!21 : In the system (3.16)-(3.17) if x a0 e B A C and p(t) is constructed as it was 
described above the states x a (t) of the system belong to B AC (i.e. x a (t)e B A)C ) for all t and for all 
r(t). 

Theorem 3.3 1121 : For the system (3.16)-(3.17) with |i(t) constructed as above the following 
is always true 

ifx a0 eB AC then lly(t)lloo ^ 1 Vinputr(t) 
if x a0 £ B a c then IlyCOlU < g(x a0 ) Vinputr(t) 

Theorem 3.4 IT 21 : At every time t^ if x a (to)e B A c , then the time-varying gain (iftg) is the 
maximum possible such gain so that 0 < Ufo) ^ °°, and Vr(t), t > to 3 p.(t), t > ^ such that the 
output lyj(t)l < 1 V i, t > to- If x a (t 0 )e B AjC then such a gain j-i(to) does not exist. 

A control structure that ensures that the control u(t) will never saturate is shown in figure 
3.3. The time-varying rate, introduced in this section, will be used as a Reference Governor (RG). 



Figure 3.3: Control structure with the RG operator 


4. Description of the Control Structure with the Operator RG 

In the proposed control structure shown in figure 3.3 the Reference Governor will mask out 
"large" references so they will not enter into the closed loop system. Choosing the Reference 
governor appropriately one can ensure that the controls never saturate so that the feedback system 
operates linearly. 

To facilitate our discussion let us assume the following models for the systems shown in 
figure 3.3. 
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Plant model: 

x(t) = Ax(t) + Bu s (t) 

(4.1) 


y(t) = Cx(t) 

(4.2) 

Compensator model: 

x c (t) = A c x c (t) + B c e(t) 

■(4.3) 


u(t) = C c x c (t) 

(4.4) 


e(t) = r^(t)-y(t) 

(4.5) 

Saturation model: 

u s (t) = sat(u(t)) 

(4.6) 


where r(t) are the reference signals, y(t) are the output signals, u(t) are the control signals 
generated by the compensator, u s (t) are the saturated (output of the saturation ) control signals. 

In addition, consider the linear closed loop system (i.e th-. system without the saturation) 
with the controls as output and assume the following representation 

x c i(t) = A cl x cl (t) + B cI r(t) (4.7) 

u (t) = C C | x cl (t) (4.8) 



X 

c 


A 

■BC' 


B 1 

C 




C 

C 



where x c i ~ 


II 

u 

< 

BC 

c 

A 

B . = 
cl 



X 


Ip 

- 


0 


c, = [ c . °] 


Following tire discussion of section 3 one can inject a time varying rate |i(t) at the inputs of a 
linear time invariant system and the outputs of that system will remain bounded. Consider the 
closed loop system (4.7)-(4.8) and assume that a time-varying rate (4.9)-(4.1 1) is introduced at the 
references as shown in figure 4.1 




Figure 4.1: The basic system for calculating |i(t). 


z(t) = p(t)e r (t) 


(4.9) 
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( 4 . 10 ) 

( 4 . 11 ) 


e r (t) = r(t) - r^(t) 
r>(0 = z(t) 

As in section 3, the time varying gain |l(t) will be chosen so that if r(t) is small enough never 
to cause control saturation then r(t) = r^(t), in contrast, if r(t) is large, then |i(t) will limit the 
references so that the controls will remain bounded. We now combine the dynamics of the rate 
limiter (4.9)-(4. 11) with the dynamics of the closed loop system (4.7)-(4.8) to obtain an 
augmented system 


x a (t) = A a x a (t)+B a |i(t)e r (t) (4.12) 

u(t) = C a x a (t) (4.13) 


where x (t) = 
a w 


z(t) " 

A = 

0 

0 

B = 

1 ‘ 

u 

0 

n 

u 

_ X clW_ 

a 

B , 

cl 

A . 

cl 

a 

0 

a 1- 


The objective here is to construct j_i(t), 0 < |i(t) < in such a way so that for any error e r (t) 
the controls u(t) never saturates. This is similar to designing a time-varying rate so that the output 
of a linear system remains bounded (section 3). At first, a function g(x) and a set B A c have to be 
defined. The symbols g(x) and B A C should be thought as generic symbols and, when they are 
used, they are always defined to avoid confusion. 



g(x a o) ; §( x ao) = ll u (0IU 

(4.14) 

where 

X a (t) = A a x a (t); x a (0) = x a0 

(4.15) 


u(t) = C a x a (t) 

(4.16) 


B A,C = { x: g(x) ^ 1 } 

(4.17) 


For the function g(x) to be finite the linear system in eq. (4.15) has to be neutrally stable. 
Even if the plant is unstable the compensator has been designed to stabilize it and the system (4. 15) 
is always neutrally stable. Consequently, the construction of |i(t) is given by: 

Construction of p(.t): 

For every time t choose |l(t) as follows 

a) if x a (t)e IntB AC then |i(t) = °° which implies that r(t) = r^(t) 

b) if x a (t)e BdB AC then choose the largest |i(t) such that 

0 < fi(t) < °° 


(4.18) 

(4.19) 
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g(x a (t)+e[A a x a (t)+B a ^(t)c r (t)]) - g(x a (t)) ^ 

lim sup — - u 

e— >o e (4.20) 

or for the points where g(x) is differentiable choose |i(t) such that 

0<p(t)<°° ( 4 - 21 ) 

Dg(x a (t))[A a x a (t)+B a |i(t)e r (t)] <0 Vt>0 (4.22) 

where Dg(x a (t)) is the Jacobian matrix of g(x a (t)) as in definition 3.2. 

c) if x a (t)£ B AjC then choose p(t), 0 < |i(t) < «= such that the expression (4.20) is 

minimum. 

The closed loop system with the RG operator has the following good properties. 

a) The controls in the closed loop system will never exceed the limits of the 
saturation and thus the direction of the control vector is not affected by the saturation. Hence, 
any inversion of the plant by the compensator is not prevented. ■ 

b) Integrators or slow dynamics in the compensator do not windup. 

The main disadvantage of this method is that the construction of |i(t) requires the 

measurement of the plant states. More research is needed to assess if estimates of the states can be 
used to approximate the real (i(t). 

It is clear that in principle this control structure can be used for any plant and any 
compensator as long as the linear closed loop system is stable (true for all sensible control 
systems). Because of the difficulty to compute |i(t), it is recommended to use the control structure 
with the operator RG in feedback system with unstable plants and/or unstable compensator only. 
For control systems with stable open loop plants and stable compensators one can use the 
methodology described in [11]. 
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r p.(t) 

I 

I 

I 



Figure 4.2: Control structure with the operator RG. 


4.1 Stability Analysis for the Control Structure With the RG 

The simple closed loop system form r(t) to r^(t) (which is part of figure 4.2) is given by the 
following: 

z(t) = |i(t)(r(t) - z(t)) (4.23) 

i>(t) = z(t) (4.24) 

The system (4.23)-(4.24) is a BIBO stable system, i.e. for bounded r(t) the signal r^(t) is 
also bounded. This can be shown formally by using Lyapunov stability theory with a Lyapunov 
function of V = z T (t)z(t) where V = 2|i(t)z T (t)(r(t) - z(t)). If r(t) is bounded the function 'V is 
negative definite for large z(t) and thus, z(t) will be bounded. 

With the RG operator the controls never saturate so the system from r^(t) to y(t) (which is 
part of figure 4.2) is a linear system; it is also assumed to be stable since one of the purposes of 
K(s) is to stabilize the linear feedback system. As a result, the control system from r(t) to y(t) is 
BIBO stable. This is an important fact because when the open loop plant is unstable the linear 
control system in the present of saturation may not be BIBO stable for all reference signals. 

Since the RG operator is outside of the closed loop system, when disturbances are present 
one cannot guarantee that the control will not saturate. In fact, there always exists a disturbance that 
will cause saturation and instability. The following analysis is done only for output disturbances 





Page 13 


d(t), similar analysis can be performed for other type of disturbances as well. 

It is clear that if r(t) is chosen so that the controls will reach the saturation limits, then there is 

a disturbance with lld(t)lloo < £, Ve > 0 such that the controls will exceed the limits of the 

x 

saturation. To avoid that, one can introduce an artificial saturation level s = [ Sj s m ] with s ; < 1 
and choose ji(t) so that the references will never cause the controls of the system to exceed the 
artificial saturation limit s. Then bounds can be defined, as we shall do in theorem 4.1, for the 
disturbances so that if the disturbances do not violate those L^, bounds the controls will always 
remain within the real saturation limits. 

In effect, the controls action can be used, partly, to track commands (lluj(t)lloo < s ; ) and, 
partly, to reject disturbances (IIUjCOIL, < 1-S;).-The artificial saturation s is "reserving" part of the 
control action for command following and the rest of the control action is used for disturbance 
rejection. In theorem 4. 1 the relationship (trade-offs) between s and the bound on the 
disturbances that can be rejected will be given. 

To ensure that only part of the control is used for command following the operator RG can be 
used to guarantee that llu^OII^ < Sj. The computation of p.(t) for this case is similar to the case 
where the saturation limit is 1. For example, in the computation of |i(t) one could scale the 
compensator so that the control saturation limits instead of s they would be 1. In the 
implementation, by rescaling the compensator, the actual saturation levels will be s again. 

Theorem 4,1 IT 21 : If the RG operator is used in any feedback system so that the controls 
(IIU|(t)lloo ^ S;), for some vector s the following is true. With zero initial conditions, the closed loop 
system with the RG operator will have bounded controls (llu(t)oo < 1) and bounded outputs for any 
reference and for output disturbances that satisfy the following condition. 

(4.25) 

where hjj is impulse response of the ij* element of the following transfer function matrix 
; H(s) = [I + K(s)G(s)]' 1 K(s) (4.26) 

From the previous discussion theorem 4.1 can be used to illustrate the trade-offs between 
"good" command following and "good" disturbance rejection. There are two ways to use theorem 
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4.1. 

(a) If we know upper bounds on the output disturbances that exist in the operating 
environment of the control system, the following is true; one can compute the artificial saturation s 
so that all possible disturbances will be rejected. These upper bounds usually come from 
experimental data and the specific operating environment of the system. Then the vector s is 
computed by the following: 

1 “ 

(4.27) 

1 

An operator RG will be included in the control system to guarantee that the references will 
never cause the controls to exceed the artificial saturation s. In this context, if S[ y for some i, is 
negative then there exists a disturbance that will cause the system to saturate even if r(t) = 0 for all 
t. If is positive, for all i, then there is a disturbance (d^(t) = ndj(t)liooSign(h|j), for some i) with 
magnitude within the specified upper bounds and some reference, to cause the controls to reach the 
limits of the real saturation (±1). In that sense theorem 4.1 is not conservative. 

(b) If the disturbances are not known then the control designer has to define the artificial 
saturation s. The value of s should be specified so that with, liu(t)lloo ^ s, there is enough control 
action for the system to perform (command following) well. One can compute s by using 
experimental data, the specifications of the control system, and the specific application. With the 
value of s one can compute an upper bound for the disturbances that will never cause saturation 
(±1) as follows: 

(4.28) 

From theorem 4.1 (expression 4.25) is evident that the smaller the disturbances (in the 
operating environment of the system) the better the command following. 

In addition to disturbances, modelling errors can cause the feedback system to saturate and 
thus degrade the performance or even to drive the system unstable. At this point, it is not clear how 
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to define a class of modelling errors so that the closed loop system with the RG operator will 
remain stable in the presence of those modelling errors. 

5.2 Simulation of the FI 6 Aircraft 


As it was described previously the introduction of the saturation in the a closed loop system 
when the open loop plant is unstable can 

(a) cause instability of the closed loop system 

(b) cause integrator windups 

(c) alter the directions of the controls and thus affect the performance of the system. 
The purpose of the next example is to illustrate these problems and to show how the new 

control design method solves these problems. 

Consider a model of the AFTI-16 (Advanced Fighter Technology Integration) aircraft, which 
is a modified F-16 aircraft. The following linear time invariant model is an approximation of the 

u s (0 

(4.29) 

(4.30) 


aircraft longitudinal dynamics at 3,000 ft altitude and .6 Mach velocity [12]. 

-2.516 -13.136 

-.1689 -.2514 

-17.251 -1.5766 

0 0 


x(t) = 


' -.0151 

-60.5651 

0 

-32.174 ' 


-.0001 

-1.3411 

.9929 

0 

x(t) + 

.00018 
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and in compact form 

x(t) = Ax(t) + Bu s (t) 
y(t) = Cx(t) 


(4.31) 

(4.32) 



u(t) 


where x(t) = 


forward velocity (ft/sec) 
a(t) angle of attack (deg) 

q(t) pitch rate (deg/sec) 

0 (t) pitch angle (deg) 


u(t) = 


y(t) = 


S c (t) 

elevator angle (deg) 

limit at 25 

5 f (0 

flaperon angle (deg) 

limit at 25° 

cc(t) 

angle of attack (deg) 


0(0 

pitch angle (deg) 



(4.33) 


(4.34) 


(4.35) 


Assume that we wish to design a closed loop system so that the F16 follows angle of attack 
and pitch attitude with zero steady state error required for step commands. Linear control theory 
will be used to design the closed loop system, then the linear design will be modified as indicated 
previously with a time varying rate |i(t). Finally, simulations will be performed to assess the 
benefits of the new design methodology. 

To obtain the linear closed loop system, integrators have to be added at the controls; and the 
augmented system is given by the following 


X aW = A a X aW +B a U aW 
y(t) = c a x a ( t ) 

u(t) = ^u a (t) 


(4.36) 

(4.37) 

(4.38) 
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A linear compensator was designed for the augmented system to control the angle of attack 
and the pitch angle tracking errors. The LQG/LTR methodology was used to design the 
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Figure 4.3: Closed loop system for the F16 example with RG. 

We now deal with a multivariable control system for an unstable open loop plant with 
integrators and a saturation element in the forward loop. Without the RG operator the control 
system is expected to have the following problems (a) for certain references r(t) the outputs y(t) 
will be unbounded, (b) integrator windups may be present, and (c) the saturation can alter the 
direction of the controls and thus degrade the transient performance of the system. 

Three types of simulation were performed for the closed loop system shown in Figure 4.3. 
These different types of simulation are the following: 

1) The first simulation is for the linear system. Again we assume that the compensator K(s) 
we designed yields desirable linear responses. 
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2) In the second simulation the saturation element is added to the linear system, u s (t) = 
sat(u(t)), without any other modifications. This simulation is referred to as the simulation for the 
system with saturation. 

3) In the third simulation u s (t) = sat(u(t)), and ja(t) is computed entirely on-line by the 
method given in section 5.2.2. The simulation was performed in a Macintosh 512K and it required 
approximately 15-16 hours. This simulation is referred to as the simulation of the system with 
saturation and the RG. 

At first, the linear system was simulated for r = [ 0 10] T corresponding to a 10° pitch angle 

with zero (trim) angle of attack. Figures 4.4 and 4.5 show the output and control responses of the 
linear system. Note that the controls have "impulsive" action at the beginning and they exceed by 
far the 25° limits so saturation is expected. Also note that the maximum flaperon control value is 
approximately 83°. 

We remark that the quality of the linear output transients (figure 4.4) are not particularly 
"nice" due to the overshoots present. However, for the sake of comparisons that follow we shall 
assume that they represent desirable shapes. 

Figures 4.6 and 4.7 show the output and control responses of the system with saturation. 
The closed loop system for the reference input r = [ 0 10] is unstable. Note that, in general, 

when the open loop system is unstable and saturation at the controls is present the resulting control 
system is only locally stable. 

Figures 4.8 and 4.9 show the output and control responses of the system with saturation and 
the RG operator. The stability of the closed loop control system is recovered and the shape of the 
transient response is similar ( but slower, as expected) to the linear response. Compare figures 4.4 
and 4.8; they are almost identical ! !!. Also, the controls u(t) never exceed the saturation limits, as 
guaranteed by the design methodology. 

Figure 4. 10 show the modified reference ty^t). Since the first reference is zero the r ml (t) is 
zero Vt and it was not plotted. Note that the r^t), commanded pitch attitude, starts at =3° and 
ramps up to the desired steady state value of 10°. The reason that the 1)^(0 is initially =3 is 
because the linear system with an r^O) of =3° will have controls with maximum at =25 
(remember that with an r m2 (t) of =10° the controls had a maximum value of =83 ). Then as the 
outputs follow the modified references the r^(t) approaches r(t) in such a way that the controls will 
never exceed the saturation limits. 


Output y(t) (deg.) 
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Output y(t) for the F 16 closed loop system vithr=[ 0 10 ] T 



0.00 2.00 4.00 6.00 8.00 10.00 

Time (sec.) 

Figure 4.4: Output response of the F16 linear system, (r = [0 10] T ). 
Control in the F 16 closed loop system vithr=[ 0 10 ] T 



Time (sec.) 

Figure 4.5: Controls in the F16 linear system, (r = [0 10] T ). 





Page 20 

T 

Output in the F 16 closed loop system vithr=[ 0 10 ] 



Time (sec.) 

T 

Figure 4.6: Output response of the F16 system with saturation, (r = [0 10] ). 

T 

Control in the F16 closed loop system vithr=[ 0 10 ] 



0.0 2.0 4.0 6.0 8.0 10.0 

Time (sec.) 

T 

Figure 4.7: Controls in the F16 system with saturation, (r = [0 10] ). 
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Output for the F 16 closed loop system vithr=[ 0 1 0] T 



0.0 2.0 4.0 6.0 8.0 10.0 

Time (sec.) 

T 

Figure 4.8: Output response for the F16 system with saturation and the RG, (r = [0 10] ). 

T 

Control for the F 16 closed loop system vithr=[ 0 10 ] 



Time (sec.) 

T 

Figure 4.9: Controls in the FI 6 system with saturation and the RG, (r = [0 10] ). 





12.0 



0.0 2.0 4.0 6.0 8.0 10.0 

Time(sec.) 


Figure 4. 10: r^ 2 (t), the commanded pitch, in the F16 system with saturation and the 
RG, (r = [0 10] T ). Insert: Blowup with 0 < t < .6. 


T 

A second simulation was performed for the same system with reference r = [2.5 2.5] . 

Now we are commanding simultaneously 2.5° angle of attack and pitch. Figures 4.1 1 and 4.12 
show the output and control responses of the linear system. Again the controls exceed the limits of 
25° and saturation is expected. 

Figures 4.13 and 4.14 show the response of the system with saturation from the output 
response one can see that the integrators windup although, now, the closed loop system is stable. 
In addition, the direction of the outputs is not similar to the direction of the outputs in the linear 
system and thus the control system does not behave as it was designed to behave. 

Figures 4.15 and 4.16 show the output and control response of the system with saturation 
and the RG. The controls never exceed the limits of the saturation and thus integrator windups are 
not present. The output response verify the absence of integrator windups. The output response is 
slower because of the limited controls but the direction of the outputs is similar to that of the linear 
system. Figure 4.17 show the modified reference r^(t). 




T 

Output for the F 16 closed loop system vith r=[ 25 .25 ] 



0.0 2.0 4.0 6.0 6.0 10.0 

Time (sec.) 

T 

Figure 4.1 1: Output response for the F16 linear system, (r = [2.5 2.5] ). 


T 

Control for the F 16 closed loop system vith r=[ 2.5 2.5 ] 



Time (sec.) 

T 

Figure 4.12: Controls in the F16 linear system, (r = [2.5 2.5] ). 
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Output for the F 16 closed loop system vith r=[ .25 .25 ] T 



0.0 2.0 4.0 6.0 8.0 10.0 

Time (sec.) 

T 

Figure 4.13: Output response for the F16 system with saturation, (r = [2.5 2.5] ). 


T 

Control for the F 16 closed loop system vithr=[ 2.5 2.5 ] 



Time (sec.) 

T 

Figure 4.14: Controls in the F16 system with saturation, (r = [2.5 2.5] ). 
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T 

Output for the F 16 closed loop system vith r=[ 2. 5 2.5] 



0.0 2.0 4.0 6.0 8.0 10.0 

Time (sec.) 

T 

Figure 4.15: Output response for the F16 system with saturation and the RG, (r = [2.5 2.5] ). 

Control for the F 16 closed loop system vith r =[ 2.5 2. 5 ] T 
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Figure 4.16: Controls in the F16 system with saturation and the RG, (r = [2.5 2.5] ). 
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5. Concluding Remarks 

In this paper it has been shown how that the operator RG can be used to design control 
systems for plants with multiple saturations. The operator RG preprocesses the reference signals in 
such a way so that the references never cause the controls in the closed loop system to saturate. 
Typically, sudden large step commands are translated by the RG operator into slower commands, 
so as to allow the limited controls not only to stabilize the system but also to eventually track the 
reference. Thus the signals in the closed loop system remain bounded for any reference and if 
integrators are present in the loop they never windup. The control structure with the operator RG 
can be used in any closed loop stable linear feedback system. 

In addition, we have shown how to define disturbance sets so that, if the disturbances belong 
to these sets, then the outputs of the system remain bounded. With this new design methodology 
one can distribute the control action among rejecting disturbances and following references as it is 
needed for specific applications. 

The main benefits of the methodology are that it leads to controllers with the following 


properties: 
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(a) The signals that the modified references never cause saturation. 

(b) Possible integrators or slow dynamics in the compensator never windup. 

(c) The closed loop system has inherent stability properties. 

These properties were demonstrated in simulations of the F8 aircraft (stable) model and an 
academic example. 
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